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Abstract 

A simple, yet accurate solution of the electron-phonon coupling problem in 
Ceo is presented. The basic idea behind it is to be found in the parametrization 
of the ground state electronic density of the system calculated making use of 
ab-initio methods, in term of sp 2+x hybridized orbitals. This parametrization 
allows for an economic determination of the deformation potential associated 
with the fullerene's normal modes. The resulting electron-phonon coupling 
constants are used to calculate Jahn-Teller effects in C^ , and multiple satel- 
lite peaks in the corresponding photoemission reaction. Theory provides an 
accurate account of the experimental findings. 
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The discovery of fullerenes |L|, and the subsequent "mass synthesis" of these hollow 
molecules, have prompted the search for a new class of materials such as fullerides (salts 
like K3C60), fullerites (Ceo molecular crystals), etc., in which fullerenes play the role of the 
building blocks. The understanding of the response of these building blocks to external fields 
is essential in the characterization of the associated new materials. A central issue in this 
quest is how the electronic properties of fullerenes are modified by the coupling of electrons 
to phonons. This question can be answered in terms of ab-initio methods (cf., e.g., Refs. 
HID and Refs. therein; cf. also Ref. ||). These methods are, however, computationally 
demanding and not particulary transparent, and much of the physical insight is lost in the 
complexity of mathematics. 

In the present paper we present a simple, yet accurate solution of the electron-phonon 
coupling problem particularly suited for fullerenes, and apply it to Cg . The central idea 
behind this solution is based on a parametrization of the ground state electronic density of 
the system calculated making use of ab-initio methods in terms of hybridized atomic orbitals. 
In particular, in the case of Cqq fullerene, of sp 2+x orbitals. Such parametrization allows for 
an economic determination of the deformation potential associated with the different normal 
modes. The resulting electron-phonon coupling constants are used to calculate Jahn-Teller 
effects in Cg arid multiple satellite peaks in the corresponding photoemission spectrum. The 
resulting cross section agrees well with the experimental data, better than any of the cross 
sections obtained making use of the electron-phonon coupling constants available in the 
literature [|3],[|,|6|-p] . Furthermore, the extension of the hybrid orbital model to describe the 
electron-phonon coupling phenomenon in fullerenes other than Cqq as well as in nanotubes 
is simple to carry out. 

The electron-phonon coupling is determined by the matrix element of the deformation 
potential 



where {R} represents the set of ionic coordinates with equilibrium values {-R }, (Qp) k the 
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displacement field associated with the normal modes of the /3-ion in the k=(x,y,z) directions, 
while V e is the total potential acting on the electrons. This potential can be expressed, in 
the local density approximation (LDA), as 

K(f, {R}) = V p L scudo (r, {R}) + V Hartrcc (p(f, {it!})) + 

vMr,{R})) + vZL(f,{R})- (2) 

The first three terms are local and arise from: a) the pseudopotential associated with the 
ions; b) the Coulomb interaction acting among the electrons (Hartree field); c) exchange 
(Fock field) and correlation effects. The last term in Eq. (|2|) reflects the non-local part of 
the ionic pseudopotential. 

Because the first and the last terms of Eq. (Q) display simply an explicit dependence on 
the ion positions, the calculation of their gradient (cf. Eq. ([]])) presents no difficulties. This 
is not the case for the second and third terms of Eq. (|2|), which depend on {R} implicitly 
via the electronic density, a fact which can be used to write 

E Qp ■ ^fiViW* {m = E ^P-Qp ■ %p(f, W)> (3) 

where i=(Hartree,xc). Because in the LDA there is an explicit relation between and 
p, the derivative dVi/dp can be calculated analytically. Consequently, the basic difficulty 
associated with the calculation of the deformation potential defined in Eq. ([!]), and thus 
of the electron-phonon coupling constants, lies in the calculation of the gradient of the 
electronic density along the normal displacements. To overcome this difficulty, we shall 
parametrize the LDA results in term of sp 2+x hybrid orbitals, in such a way that three of 
the four orbitals are directed along the carbon bonds of fullerene Cgo, while the fourth takes 
care of the 7r-bonding present in the hexagons and is directed essentially perpendicularly to 
the fullerene surface. 

In keeping with the fact that in Cgo each atom has three nearest neighbours, the hybrid 
orbitals we are interested in can be written as 

|0i > = ai|s > +h\p x > +ci\p y > +d 1 \p z >, 
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|02 > = a 2 \s > +5 2 {m 2 \p x > +n 2 \p y > +t 2 \p z >), (4) 
|03 > = a 3 \s > +5 3 (m 3 \p x > +n 3 \p y > +t 3 \p z >), 
|04 > = a A \s > +5 A (m A \p x > +n A \p y > +U\p z >). 

Here rrij = cosaj, rij = cosPj, tj = cos-yj, and aj, j3j, jj (j=2,3,4) are the angles which 
define the direction of the bond between a Carbon atom and each of its three nearest 
neighbours, in a system of reference centered on the atom. After having fixed the direction 
of these orbitals, there still remain ten free parameters in Eq. parameters which are 
completely determined by the orthonormalization condition. 

To describe the radial dependence of the | s) and | p) orbitals, we have used the functions 
R s = -^=e^~ r l a ^ and R v = r — re <y ~ r / (72 \ usually employed in the description of the Carbon 

\/<T? F \/3<tS 



atom. We have however adjusted the parameters 0\ and a 2 in order to obtain the best fit 
to the LDA Ceo-density (see, for instance, Ref. 0). One can then write the contribution to 
the total density arising from a single atom and, with the help of standard techniques, carry 
out a multipole expansion of this contribution around the center of the molecule. Adding 
the contributions of the 240 electrons one obtains the total density. 

In Fig. [I] we display the two lowest multipoles of the C 60 ground state density, calcu- 
lated in the LDA including exchange-correlation effects according to the parametrization of 
Perdew and Zunger IfLOfl . The role of the Carbon atoms were taken into account in terms of 



norm-conserving pseudopotentials 0] . In the same figure we show the results of the hybrid 



orbital model, for a x = 0.78 A and a 2 = 0.31 A JLJJ. 

The next step consisits in the calculation of the gradient of the density and, through Eqs. 
(0) and (f|), the deformation potential. Within the hybrid orbital model, moving around the 
ions change the direction of the orbitals but not their shape. This means that the weights 
Sj in Eq. (^) are fixed, and the only quantities which change are the angles <x,-, (3j and 7j. 
The calculation of Vp(r, {R}) becomes then quite simple. 

In Fig. |2| we display the two lowest, and most important, multipole contributions of 
the local part of the deformation potential of Ceo for the lowest A g phonon, corresponding 



essentially to a breathing mode of the system, as calculated in the LDA. The wavenumber 
of this mode is equal to 491 cm -1 and to 496 cm -1 for Cqq in solution and in the solid 
phase, respectively [|13|,[14]]. The wavenumber and zero point motion of the isolated molecule, 
calculated in the bond-charge model, are 493 cm -1 and 53.4 10 _3 A, respectively (|IIJ. As 



seen from Fig. ^|, the hybrid orbital model provides an overall account of the ab-initio results. 
The discrepancies observed for small values of r for the multipole (L,M) = (0,0) have little 
influence on the corresponding matrix elements, because the electronic wavefunctions are 
quite small around the origin of the molecule. 

Making use of these results, one can calculate the electron-phonon coupling matrix el- 
ements in Ceo- These are the matrix elements needed, for example, in the evaluation of 
the lineshape of allowed as well as of forbidden electronic transitions in Ceo- To carry out 
similar calculations in Cg one should employ a deformation potential which is evaluated 
by making use of the electronic density of the negative ion, at the ground state geometry. 
Because the density of the 240 valence electrons of C 60 is not appreciably altered by adding 
one more electron, one expects the deformation potentials associated with Cqo and Cg~ to 
be quite similar. In fact, we have carried out fully relaxed, ab-initio calculations of the 
matrix elements in Cg and found that they agree with those of Ceo within less than 10%. 
In keeping with these results, the electron-phonon matrix elements calculated starting from 
the Ceo elecronic density and ground state ionic configuration are used in the following. 

In Ce the state t lu is occupied with a single electron. This level is separated by an energy 
of the order of the eV from neighbouring levels, while the electron-phonon coupling matrix 
elements to be found below are of the order of the meV. Consequently, it seems justified to 
consider, within the present context, that the electronic motion is confined to the ti u level. 
Under such circumstances, and because of symmetry reasons, the only possible couplings 
are to phonons with A g and H 9 symmetries (cf., e.g., Ref. |7|]). 

The matrix elements {ti u u \V def \ ti u ), where v stands for the quantum numbers of the 
phonon, are related to the partial electron-phonon coupling constant \ U /N(0) = agl/u>„, 
according to (t\ u v \V det \ ti u ) = (g v /2) W l nm (cf., e.g., Ref. |16| and Refs. therein). In the 



above expression N(0) is the density of levels at the Fermi energy, a is equal to 1/3 for A g 
phonons and 5/3 for H 9 phonons, while uj u is the energy of the phonon. The quantities W l nm 
are geometric coefficients, the index I distinguishing between the different degenerate states 
of each phonon (H 5 is five-fold degenerate while A g is single-fold degenerate). 

In Table | we display the multipole contributions to the matrix element {t\ u A g \ V de{ \ ti u ), 
associated with the lowest energy A g mode calculated making use of the LDA and of the 
hybrid orbital model. The different contributions of the model display the same sign and 
similar order of magnitude as those calculated in LDA, while the summed contribution agree 
within 20%. 

Following the same steps as those leading to the results displayed in table d the different 
matrix elements (ti u u \V de{ \ t\ u ) (z/=A g ,H 5 ) have been calculated. Our results have been 
compared with those from other theoretical calculation available in the literature [||,|],[7|||, 
as well as with the empirical values obtained from Gunnarsson's systematic analysis of 
the photoemission spectra of C^ [[J. While this analysis indicates that the coupling of 
the t iM electron to the H 5 (2) leads to the largest value of \ U /N(0), the results reported 
in Refs. []3]||{(]>@] indicate the coupling to the H 9 (7)-phonon to be the most important. In 
the hybrid orbital model discussed above, the largest coupling of the t iM level is to the 
H 3 (2)-phonon, in agreement with the analysis by Gunnarsson and co-workers [[Hf . 

Making use of the matrix elements (ti u v \V dcf \ ti u ) (z/=A g ,H g ) calculated in the hybrid 
orbital model, and of the results of the bond charge model [[n| to describe the prop- 



erties of the phonons, we have solved the total electron-phonon Hamiltonian contain- 
ing an electronic term, a phonon term, and a linear coupling term (cf., e.g., Ref. ||) 
in a basis of one ti u electron and up to three phonon states. The lowest eigenvalue 
I tiu) — C | ti u ) + J2u Cy 1 \t\u® v ) + • • • was calculated using the Lanczos method. The 
first term in describes a state with no phonons, the second term a state with one 

phonon, etc. Making use of these results, we have calculated the photoemission cross sec- 
tion assuming the sudden approximation, where the emitted electron does not interact with 
the system left behind ||. 
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In Fig. |3| we show the results of the hybrid orbital model, in comparison with the re- 
sults of the analysis of the photoemission data carried out by Gunnarsson and co-workers 
||. Although the hybrid orbital model leads to a somewhat weaker electron-phonon cou- 
pling than required by the experimental finding, and consequently to a somewhat too large 
value of C(°\ it provides a much better account of the empirical spectrum than the other 
theoretical models, whose partial electron-phonon coupling constant have been reported in 
Refs. 1,1,01 (cf. also table I of Ref. [g). 

We conclude that the hybrid orbital model of the electron-phonon coupling displays a 
number of attractive features: i) it leads to matrix elements of the deformation potential 
which reproduce quite accurately the results of ab-initio calculations; ii) it provides an 
excellent account of the photoemission spectra of C^ ; iii) it is quite economic to use, and 
can be extended at profit to fullerenes other than Cqq as well as to nanotubes. 

Finantial support by NATO under grant CRG 940231 is gratefully acknowledged. We 
are also indebted for support to INFM Advanced Research Project "CLASS" . 
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TABLES 

TABLE I. Different multipole contributions to the matrix element {t\ u A g \ V dci \ ti u ), associated 
with the lowest energy A g mode of Ceo {tiu is the LUMO of the molecule). 

Matrix Element [meV] 



L= 


LDA 


Hybrid Orbital Model 





-47.739 


-49.506 


6 


-3.319 


-2.369 


10 


19.741 


20.275 


12 


-1.030 


-2.828 


16 


6.082 


5.869 


18 


14.940 


15.473 


20 


2.139 


1.777 


Total 


-9.399 


-11.309 



10 



FIGURES 

FIG. 1. Comparison between the results of LDA (full lines) and Hybrid Orbital Model (dashed 
lines) for the two main multipole contributions to the ground state density of Ceo- (a) and (b) 
refer to the (L,M)=(0,0) and (6,0) contributions respectively. 

FIG. 2. Same as Fig. [I] for the two main multipole contributions to the deformation potential 
associated with the lowest A g phonon of Cgo- 

FIG. 3. Results for the photoemission spectrum of Cg , obtained making use of the elec- 
tron-phonon matrix elements calculated in the present paper as well as in previous theoretical 
works The solid curve correspond to the experimental results Q. 
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